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ABSTRACT 

We have observed a sample of 14 nearby {z ~ 0.03) star-forming blue com¬ 
pact galaxies in the rest-frame far-UV ( 1150 — 2200 A) using the Cosmic 

Origins Spectrograph on the Hubble Space Telescope. We have also generated 
a grid of stellar population synthesis models using the Starburst99 evolutionary 
synthesis code, allowing us to compare observations and theoretical predictions 
for the SiIV_1400 and CIV_1550 UV indices; both are comprised of a blend of 
stellar wind and interstellar lines and have been proposed as metallicity diag¬ 
nostics in the UV. Our models and observations both demonstrate that there 
is a positive linear correlation with metallicity for both indices, and we find 
generally good agreement between our observations and the predictions of the 
Starburst99 models (with the models slightly under-estimating the value of the 
indices due to contributions from interstellar lines not simulated by a stellar 
population synthesis code). By combining the rest-frame UV observations with 
pre-existing rest-frame optical spectrophotometry of our blue compact galaxy 
sample, we also directly compare the predictions of metallicity and extinction 
diagnostics across both wavelength regimes. This comparison reveals a corre¬ 
lation between the UV absorption and optical strong-line diagnostics, offering 
the first means of directly comparing ISM properties determined across different 
rest-frame regimes. Einally, using our Starburst99 model grid we determine the¬ 
oretical values for the short-wavelength UV continuum slope, /Sig, that can be 
used for determining extinction in rest-frame UV spectra of star-forming galax¬ 
ies. We consider the implications of these results and discuss future work aimed 
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at parameterizing these and other environmental diagnostics in the UV (a suite 
of diagnostics that could offer particular utility in the study of star-forming 
galaxies at high redshift) as well as the development of robust comparisons 
between ISM diagnostics across a broad wavelength baseline. 


1. Introduction 


A galaxy’s redshift plays a dominant role in determining the rest-frame wavelength 
bands in which it is most easily observed. The rest-frame UV spectra of local galaxies 
can be observed with UV-sensitive instruments. The rest-frame UV spectra of galaxies at 
z ~ 2 — 3, on the other hand, must be observed in the optical. Such optical observations 
offer many practical advantages over observing the more distant galaxies in the IR, in 
order to access their rest-frame optical spectra, due in large part to an abundance of 
strong night sky emission lines in the IR. With z ^ ^ ~ 3 galaxies sampling the regime 
when star formation density peaks (iHopkins fc BeacomI 120061 ). it is therefore advantageous 
to utilize UV data from local starburst galaxies to develop spectral interstellar medium 
(ISM) diagnostics that can ultimately be applied to optical observations of higher-redshift 
galaxies. 

Historically, the development of diagnostics for determining oxygen-abundance metal- 
licity (log(0/H) -|- 12) in star-forming galaxies has largely been confined to the optical. 
The “direct” metallicity diagnostic uses the auroral [OIII] A4363 line to determine the 
electron te mperature, Tp, of the gas in galaxy, which can then be converted to a metal¬ 
licity (e.g. lAlleiill984l ). However, this auroral line is weak even at low metallicities and 
this method is therefore rarely usable in high-metallicity or faint galaxies. As a result, 
diagnostic s utilizing observed st rong-line ratios have been calibrated against Tg metallic¬ 
ities (e.g. iPettini &: Pagell l2004lf . Strong-line diagnostics have also been developed based 
on theoretic al calibrations using populat ion synthesis models along with photoionization 
models (e.g. iKobulnickv &: Kewlevll2004^. or th rough a combination of observational and 
theoretical methods (e.g. iDenicolo et al.M2002r). Unfortunately, the numerous metallic - 
ity diagnostics available in the optical (jKobulnickv &: Kewlevll2004 iPettini &: Pagelll2004l l 
have been found to show poor agreeme nt, although conversions do exist which facilitate 
comparisons between different methods (jKewlev fc EllisonI 120081 ). 


Fewer diagnostics are available in the UV. Interstellar absorption lines in the UV 
are produced by the resonance transitions of abundant ions. These lines are stronger in 
starbursts, due to the larger interstellar velocity dispersions. The stronger interstellar lines 
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tend to be optically thick and t herefore their equivalent width s depend more on velocity 
dispersion than column density (iGonzalez Delgado et al.lll998l ). Because the more metal- 
rich galaxies tend to host more powerful starbursts, those galaxies tend to have more energy 
in the ISM. This inc reases the velocity di spersion and therefore the equivalent widths of 
th e interstellar lin e s (ILeitherer et al.l 120111!. This m etallicity dependency has been shown 
bv iHeckman et al.l (119981 ) and ILeitherer et al.l (|2nill ). Stellar photospheric absorption lines 
can also be used as metallicity diagnostics when the spectra have the hig h signal-to-noise 
and m oderate resolution required to me asure equivalent widths of ~ 1 A ( Leitherer et al. 


201in. To aid in this r esolution problem, iRix et al.l (|2004l l employed the indices defined by 
Leitherer et al.l ( 20011 ) to determine that a region between 1935 A and 2020 A, containing 
many Fe III photospheric lines, scales with metallicity. 

Stellar w ind lines — resonanc e lines produced in the radiatively driven winds of O 
and B stars (jLeitherer et al.l l201ll ) — offer the most robust means of directly tracing 
metallicity in the UV. The lines originate mainly from gas at T > lO^K which has been 
i onized by radiation from massive stars along with collisional processes in the outflow 
([ShaDlev et al.ll2003l ). The winds are driven by photon momentum t ransfer through meta l¬ 
line absorption , so th e stellar wind lines are metallicity-dependent (ILeitherer et al.ll2011h. 


1984 


Faber et al 


Leitherer et al.l (l201ll) d eveloped a set of indices, similar to the Lick indices (jBurstein et al. 


19851 ). for the purpose of measuring interstellar and stellar wind lines 


in the UV. These indices measure the equivalent width of blended absorption lines in a 
systematic manner. This allows for the consistent measuring of line strengths even in 
spectra where the resolution is not high enough to measure individual line widths. 

Unfortunately, robust conversions between UV and optical metallicity diagnostics are 
non-existent. The optical and UV light of star-forming galaxies originate from different 
processes — UV stellar spectra consist mainly of the light from O and B stars of young 


stella r populations, while optical stellar spectra trace older stellar populations (jRix et al 


20n4h — making a clear relationship between diagnostics in the two regimes even more 
critical. 

A similar difficulty exists for spectral diagnostics of dust extinction in star-forming 
regions. The effective extinction is defined as th e difference in magnitudes of the attenuated 
light and the dust-free source (jCalzettil 120011 ) , and sensitive diagnostics are required to 
determining the intrinsic attenuation due to dust. In the optical, the Balmer decrement 
(Ha/H/3) is most commonly used to quantify extinction effects. The intrinsic Balmer 
decrement is expected to be 2.86 assuming case B recombination (Tg = 10^ K and rie ~ 


10^ — 10^ cm following IOsterbrocklll989l ). with a higher observed decrement indicating 


increased dust content. When UV and optical spectra are both available, line-to-continuum 
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ratios Ha/UV and H/5/UV are also popular (jCalzettill200lh . 


The slope of the UV continuum is another widely-used extinction diagnostic. The UV 
spectral slope, /? is determined by htting the continuum to the function /(A) oc A^, where 
/(A) is the flux density, in a specified wavelen gth range. Th e wavelength range 1250~2600A 
(which defines j32&) is most commonly used (ICalzettill200ll) . The UV spectral slope traces 
the reddening effects of dust on the continuum. Dust is proportional to metallicity and 
gas content, which evolve in opposite directions with time. Whereas metallicity increases 
with time as more stars have formed and produced metals, more gas is locked away in stars 
as time goes by. Thus the dust column density peaks at som e point within the galaxy ’s 
evolution, with an overall peak predicted at a redshift of ^ ~ 2 (ICalzetti fc Heckmanlll999l ). 


However, there are degeneracies between extinction effects and other galaxy properties 
including the age of the stellar population, which dominates the spectral energy distribu¬ 
tion (SED); the effects of metallicity; and the initial mass function (IMF) . These become 


espec ially troublesome when only low-resolution spectroscopy is available (jPapovich et al. 


20011 ). Furthermore, reddening effects are sensitive to the geometry of the dust distribution 
when dealing with extended sources such as star-forming regions. Ty pically the light is not 
as attenuated as would be expected from a simple dust-screen model (IWitt et al.lll992l ). In 
addition, spectral coverage in the UV often spans a narrower wavelength range than that 
defined by the /326 diagnostic, limiting its efficacy as an extinction diagnostic at higher 
redshifts. 


In this paper we present a preliminary analysis of UV ISM diagnostics for metallicity 
and extinction. Using UV and optical spectra of nearby {z ~ 0.03) star-forming galaxies 
(Section 2) along with a grid of stellar population synthesis models (Section 3), we examine 
both observed and theoretical predictions for the utility of two potential UV diagnostics for 
metallicity and their correlations with optical strong-line diagnostics (Section 4). We also 
present our measurements of the shorter-wavelength UV spectral slope f3is and discuss its 
utility as an extinction diagnostic in both nearby and high-redshift star-forming galaxies 
(Section 5). Finally, we discuss our conclusions and plans for future work (Section 6). 


Observations 


We originally selected a sample of 31 nearby (0.003 < 2 : < 0.029) star-forming blue 
co mpact galaxies ( BCG s) wit h optical spectra from a larger set of 97 galaxies presented 


m 


Kong fc Chend (j2002l ) and iKong et al.l (j2002l ) . Rest-frame optical spectra of the galax¬ 


ies’ bright nuclei were obtained with the OMR spectrograph on the 2.16m telescope of 
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the Beijing Astronomical Observatory, using a 2""3''slit width that was a diusted based on 
seeing conditions and, where possible, aligned with the parallactic angle. iLevesque et ah 


(j2010l i then selected 36 of these BCGs (res tricting the sample to star-forming galaxies 
based on the criteria of iKewlev et al.l l2006l l , and used optical emission line diagnostics 
to determine a number of IS M properties. These in c lude two different metallicity diag- 
nostics: the calib ration of iKobulnickv fc KewlevI ([200J) and the 03N2 calibration of 


Pettini k. Pagell (120041'): for mor e discussion see Section 4.1. Subsequent study revealed that 


5 of the iLevesaue et al.l ([20101) BCGs are actually pairs of galaxies; these were excluded 
from our final sample. 


17 of these star-forming galaxies have been successfully observed by the Cosmic Origins 
Spectrograph (COS) onboard the Hubble Space Telescope (HST) as part of SNAP program 
^13481 (PI: Levesque). We used the low-resolution G140L grating at the llOSA central 
wavelength setting to provide continuous spectral coverage over the range 1150-2200 A 
at a spectral resolution of i? ~ 1800 {Av ~ 200km s“^). Each observation consisted of four 
exposures at different focal plane (FP) positions settings to dither instrumental features 
over the detector and increase overall signal to noise. Total exposure time per target varied 
from less than a minute to several kiloseconds; exposure times for each observed galaxy 
in our sample are given in Table 1. As of this writing the completion rate of this SNAP 
program is 55%; one galaxy (2MASX J05554264-I-0323322) was observed with a truncated 
exposure time of 87 seconds due to drift from single star guiding effects, and two additional 
galaxies (III Zw 43 and Haro 3) also had insufficient S/N for our analyses; these data were 
not included in this work. Our final sample of 14 galaxies, including COS exposure times, 
is given in Table 1. 


The exposures were obtained from the Mikulski Archive for Space Telescopes (MAST). 
The exposures were binned by 3 pixels (~ 1/2 of a resolution element) to improve S/N 
prior to cross-correlation and coaddition. The calibrated, one-dimensional spectra for each 
tar get were next al i gned and coadded using the custom IDL procedures described in detail 
by 


Danforth et al. 


(201C 


/]. Briefly, each exposure of a given target was aligned by cross- 


correlating on a strong interstellar absorption feature. The data were then interpolated to 
a common wavelength scale and the flux and error at each wavelength were combined using 
an exposure-time weighted algorithm. Pixels with instrumental features were de-weighted 
or excluded from the final coaddition. For our analyses the reduced spectra were adjusted 
to rest-frame wavelengths (redshifts are given in Table 1). 


^IDL routines available at http://casa.Colorado.edu/$\sim$daiiforth/costools.html 
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Fig. 1.— HST COS spectra of our star-forming galaxies. Spectra have been corrected to 
rest-frame; the strong emission features near <1215A and <1302A represent blueshifted 
geocoronal Lya and Oi, respectively. Red regions of the spectra illustrate the wavelengths 
of the Leitherer SiIV_1400 and CIV_1550 diagnostic indices; see Section 4.2 for further 
discussion. The spectra are plotted in log space to best illustrate our fits (blue line) of 
the UV spectral slope, Pis (where f(A) oc for 1250 < A < 1800 A; see Section 5). 
Variations as determined from adding and subtracting the error spectra are shown in gray. 
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The fully-combined data vary wildly in terms of data quality due mostly to source flux 
and exposure time. Since many of the star-forming galaxies in our sample are extended 
objects which fill the COS aperture, resolution for some targets may be degraded from the 
nominal R ~ 1800 typical of point sources observed with COS/G140L; however, BCGs are 
dominated by their central clusters and thus appear close to point-like, yielding only min¬ 
imal loss of resolution in COS observations. It should also be noted that the 2.5'Miameter 
field s top of the COS Pr imary Science Aperture is comparable to the 2''~3''slit widths used 
in the lKong et al.l (j2002l i observations; by centering COS on the bright nuclei of the BCGs 
we sample c ompa rable galaxy regions with our FUV spectra to those observed with the 


Kong et al.l (|2002l l optical spectra. Our spectra are shown in Figure 1. 


3. Stellar Population Synthesis Models 


Starburst99 (ILeitherer et al.lll999l : IVazauez fc Leithereiil2005l : iLeitherer &: Chenll2009l : 
Leitherer et al.ll2ni4l ) is a stellar population synthesis code which produces model spectral 
energy distributions as a function of metallicity, star formation, initial mass function (IMF), 
age, and evolution of the simulated stellar populations. This is done by combining model 
stellar atmospheres and spectra with grids of stellar evolutionary tracks. 

For this work we used this code to generate a grid of model ionizing spectra with 
varied stellar evolutionary models, metallicities, IMFs, and treatments of the star formation 
history. We used three differen t sets of evolution tracks produced by the Geneva group. 
The tracks from iMevnet et al.l (j 19941 ) (hereafter M94) simulate the evolution of a non¬ 
rotating stellar popula tion with enhanced m ass-l oss rates. We a lso ad opt the rotating and 
non-rotating tracks of lEkstrom et al.l (120121 1 and lCeorgy et al.l (j2013l h the rotating tracks 
(hereafter ROT), include the effects of stellar rotation with initial rotation rates of u = 0.4 
of the break-up ve locity, which is the peak of the velocity distribution in young B stars 
()Huang et al.l 120101 1. while the non-rotating tracks (hereafter NOROT) simulate identical 
initial conditions but do not include rotation. 


The metallicities of our models were based on the assumed abundances for the evolu¬ 
tion tracks. The M94 tracks span hve fractional heavy element abundances by mass {Z = 
0.001, 0.004, 0.008, 0.020, and 0.04), while the ROT and NOROT tracks are available at 
Z = 0.002 and Z = 0.014. The abundance for the optical high resolution spectrum in 
Starburst99 is set separately, and was chosen for each simulation to most closely match the 
abundance of the chosen stellar evolutionary tracks. 


We initially adopted a Kroupa IMF (jKroupal l200ll ). a power-law in dN/dM with 














































a = 1.3 for the 0.1 — 0.5 Mq mass range, and a = 2.3 for the 0.5 — 100 Mq mass range. 
In addition, top-heavy and flattened IMF models were made with high-mass slopes set to 
a = 1.3 and a = 3.3 while keeping the low-mass slope set at a = 1.3. As massive O and 
B stars are the main contributors to the UV spectrum, we are primarily concerned with 
the high-mass slope. Unless otherwise stated, all models presented here assume a Kroupa 
IMF. 


Our models simulated both an instantaneous burst of star formation and a region of 
continuous star formation. In models with an instantaneous burst of star formation, a 
zero-age stellar population of 10® Mq is simulated, allowing us to trace the effects of a 
single coeval stellar population; we modeled the evolution of the stellar population up to 
20 Myr in 0.5 Myr increments. The continuous star formation models adopt a constant 
rate of 1 Mq yr“^, representing the other extreme of star formation history. Here, while 
we also model the evolution of the population in 0.5 Myr increments, we consider only the 
models at an age of 5 Myr, the age at which a continuously star-forming stellar population 
reaches equilibrium (at younger ages this equilibrium is not yet reached; at older ages 
the re is little to not ev olution of the stellar FUV ionizing spectrum. For more discussion 
see 


Kewlev et al.l 1200111 . 


4. Metallicity Diagnostics 


4.1. Optical 


Using the existing optical spectra for our observed galaxies from lKong fc Chend ([2002) 
and iKong et al.l ([2002) we can calculate the metallicities (defined here as the oxygen abun¬ 
dance 12 -|- log(0/H)) using several common strong-line optical diagnostics. 


Kobulnickv fc Kewlev! (l2004ll fher e after KK04) used the stellar evolution and pho¬ 
toionization grids of [Kewlev fc Dopital (|2002 ) to derive their metallicity diagnostic. Like 


many optical metallicity diagnostics, this method depends on the ([O ii] A3727 -|- [O iii] 
AA4959, 5007)/H/3 line ratio, known as i? 23 - Unfortunately R 23 is double valued with 
metallicity; KK04 uses the additional line ratio [N ii] A6583/[0 ii] A3727 to break this 
degeneracy and determine whether to use the upper or lower metallicity result. The KK04 
diagnostic works by iterating between calculating ionization parameter and metallicity un¬ 
til the values converge. The ionization parameter formula depends on metallicity and the 
[O III] AA4959, 5007/[0 ii] A3727 line ratio. The metallicity formulae depend on the ion¬ 
ization parameter and R 23 , and have upper and lower branch versions. If the lower branch 
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gives a higher metallicity than the upper branch, the metallicity is set to 12 + log(0/H) 
= 8.4, the cutoff between the two branches. The estimated accuracy of this diagnostic is 
~ 0.15 dex. 


Pettini k. Pagell (j2004l i developed two optical metallicity diagnostics using Tg-based 


metallicities and observations of H ii regions. The first of these diagnostics (hereafter 03N2) 
is a linear fit which depends on the line ratio log[(0 iii] A5007/H/3)/([N ii] A6583/Ha)]. It is 
valid for log[(0 iii]/H/3)/([N ii]/Ha)] < 2. The estimated accuracy of the 03N2 diagnostic 
is ~ 0.14 dex. The second diagnostic (hereafter N2) is a cubic fit which depends on the line 
ratio [N ii] A6583/Ha. It is valid for —2.5 < [N ii]/Ha < —0.3. The estimated accuracy 
of the N2 diagnostic is ~ 0 .18 dex. For a more det ailed discussion of optical metallicity 
diagnostic development, see lKewlev fc EllisonI (j2008l ). 


4.2. UV - Leitherer Indices 


Leitherer et al.l (j201l|) established a set of UV indices for star-forming galaxy spectra. 
Each index is a blend of lines, mostly of interstellar or stellar wind origin. Three wavelength 
regions are defined for each index: a blue continuum, a central index bandpass, and a 
red continuum. Within each continuum range the median flux in calculated, as well as 
the midpoint in the wavelength range. The continuum is defined as the line connecting 
the representative blue continuum and red continuum points. The equivalent width is 
calculated as the sum of the difference between the data and the continuum in each bin 
within the central bandpass, normalized by the continuum. The error was determined by 
calculating the standard deviation of the blue and red continua data from the continuum 
line. This value was multiplied by the square root of the number of bins in the central 
index bandpass to find the error of the equivalent width of the index. 

Of Leitherer’s indices, SiIV_1400 and CIV_1550 are the only two dominated by stellar 
wind lines. The main feature included in SiIV_1400 is the doublet Si iv A1393,1402, while 
the main feature included in CIV_1550 is the doublet C iv A1550,1548. Both the Si iv 
A1400 doublet and the C iv A1550 doublet have interstellar components, which sometimes 
provide a significa nt contribution to the li ne equivalent wi d ths, i n addition to their stellar 
wind components (ILeitherer et al.ll201lh . iHeckman et al.l (119981 1 finds a correlation with 
dust extinction and reddening of the vacuum-UV continuum for these high-ionization lines, 
but a weaker correlation that those found for low-ionization lines. 


In their spectra. ILeitherer et ahl (j201l|) found that, by visual inspection, the interstellar 
component in the Si iv doublet could account for ~ 2/3 of the equivalent width, whereas 
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the interstellar component in the C iv doublet could account for ~ 1/3 of the equivalent 
width. While the stellar wind co mponent in the Si i v doublet only becomes apparent for 
blue giants and supergiant stars (IShaplev et al.ll2003ll. wind e ffects can be seen in the C iv 
doublet of early-0 main-sequence stars (jLeitherer et al.ll200ll) . 


4.3. Comparison of Optical and UV Diagnostics 


We measured the SiIV_1400 and CIV_1550 indices in all of our HST COS spectra (see 
Table I4.3p ; the wavelength ranges of the central bandpasses for both indices in each of our 
galaxies are illustrated in Figure 1. For comparison purposes the COS data were binned by 
a factor of 2 to approximate the O.dA resolut ion of the St a rburs t99 spectral output. We also 
used the optical emission line fluxes given in iKong et al.l (j2002l l to determine optical strong 
line metallicities using the KK04, 03N2, a nd N2 diagii o stics for each of our galaxies based 
on the optical emissio n line fluxes given in iKong et al.l (120021 ) and the analyses detailed in 
Levesque et al.l ()2010l ). 


In Figure 2 we compare the two Leitherer indices to the three optical diagnostic 
metallicities. There are statistically robust positive linear relations (Spearman’s p > 0.6, 
p < 0.02) between the SiIV_1400 and CIV_1550 indices and all three of the optical diag¬ 
nostic metallicities. 

We measured the SiIV_1400 and CIV_1550 Leitherer indices in our full grid of Star- 
burstOD models using the same definitions and calculation procedures detailed in Section 
4.2. The indices were measured from the non-normalized .ifaspec Starburst99 output 
spectra; these are theoretical UV spectra calculated between 900 -3000A at a resolution 
of ~0.4A, adopting WM-Basic models for the O star at mospheres (Pauldrach et al.ll200in 
and the Potsdam mo dels for Wolf-Rayet at mospheres (jHamann &: Grafened l2004l L for a 
detailed discussion see iLeitherer et al.l (j201fll ) . 


It should be noted that, as an evolutionary synthesis code, Starburst99 alone does not 
include any treatment of photoionization or contributions from the ISM. As a result, the 
interstellar contributions to the Leitherer indices seen in observations are not present; the 
models only simulate the stellar components (wind and photospheric) of these indices. 


The index definitions of ILeitherer et al.l (j201ll l also include an inherent metallicity 
bias. The definitions of the indices’ continuua are themselves metallicity dependent and 
thus measure a “continuum” level that is lower at high metallicities due to line blanketing 
effects, leading to a slight decrease in the equivalent widths of indices measured at high 
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Table 1. Blue Compact Galaxies Observed with HST 


Galaxy 


RA 

Dec 

z 

B 

log(0/H)+12 
KK04 03N2 

a 

N2 

SFR^ 
(Mo yr-i) 

Exp. Time 

(s) 

III Zw 12 

00 

47 

56.3 

+22 22 22.5 

0.019 

15.0 

8.57 

8.52 

8.51 

1.58 

720 

Mrk-5 

06 

42 

15.9 

+75 37 40.0 

0.003 

17.0 

8.64 

8.11 

8.13 

0.03 

1080 

Mrk-390 

08 

35 

33.1 

+30 32 03.2 

0.025 

15.0 

8.41 

8.42 

8.42 

9.31 

1568 

Mrk-105 

09 

20 

26.3 

+71 24 16.0 

0.013 

16.0 

8.94 

8.56 

8.58 

0.86 

2060 

Haro 2 

10 

32 

32.0 

+54 24 03.5 

0.005 

13.7 

8.72 

8.41 

8.39 

1.03 

420 

Haro 25 

10 

48 

44.2 

+26 03 13.2 

0.026 

15.0 

8.76 

8.33 

8.36 

11.08 

568 

Haro 4 

11 

04 

58.5 

+29 08 22.1 

0.003 

15.5 

8.22 

7.97 

7.96 

0.02 

260 

I Zw 53 

13 

13 

57.7 

+35 18 52.8 

0.016 

16.8 

8.66 

8.55 

8.62 

0.56 

2200 

H Zw 70 

14 

50 

56.5 

+35 34 18.2 

0.004 

14.5 

8.09 

8.14 

8.09 

0.08 

80 

I Zw 159 

16 

35 

21.1 

+52 12 52.8 

0.010 

15.6 

8.56 

8.25 

8.26 

0.52 

620 

I Zw 191 

17 

40 

24.8 

+47 43 58.6 

0.019 

14.8 

8.92 

8.73 

8.72 

3.23 

1120 

IV Zw 93 

22 

16 

07.7 

+22 56 32.7 

0.012 

15.7 

8.40 

8.26 

8.22 

0.57 

1120 

Zw 2220 

22 

23 

00.6 

+30 55 23.0 

0.023 

14.0 

8.88 

8.59 

8.71 

4.90 

380 

IV Zw 149 

23 

27 

41.1 

+23 35 20.2 

0.011 

13.2 

8.47 

8.41 

8.49 

0.73 

140 


^Based on analyses in lLevesgue et ^ ll2010f) . adopting 


raw emission 


line fluxes from iKong et all 120021) and 


correcting for the total line-of-sight extinction. 


Table 2. Measured Spectroscopic Properties of BCGs Observed with HST 


Galaxy 

SiIV.1400 

(A) 

C IV.1550 

(A) 

/3i8' 

a 

T 

l 

B 

a,b 

HI Zw 12 

5.41 

± 

2.74 

8.56 

± 

2.98 

-0.12 

+ 

0.14 

0.25 

+ 

0.04 

Mrk 5 

1.47 

± 

1.64 

4.97 

± 

2.10 

-1.27 

+ 

0.06 

-0.04 il 

I 0.01 

Mrk 390 

3.60 

± 

0.77 

5.83 

± 

1.14 

-0.92 

+ 

0.05 

0.47 

+ 

0.04 

Mrk 105 

4.74 

± 

0.93 

9.31 

± 

1.10 

-1.26 

+ 

0.04 

0.42 

+ 

0.04 

Haro 2 

5.11 

± 

0.48 

7.98 

± 

0.67 

-1.72 

+ 

0.03 

0.43 

+ 

0.03 

Haro 25 

3.39 

± 

2.20 

4.74 

± 

2.24 

0.87 + 0.13 

0.26 

+ 

0.03 

Haro 4 

1.45 

± 

1.37 

3.00 

± 

1.78 

-2.01 

+ 

0.08 

0.01 

+ 

0.01 

I Zw 53 

5.48 

± 

1.23 

6.65 

± 

1.42 

-1.13 

+ 

0.08 

0.65 

+ 

0.18 

H Zw 70 

4.47 

± 

1.31 

2.70 

± 

1.98 

-1.63 

+ 

0.10 

-0.08 il 

I 0.03 

I Zw 159 

2.34 

± 

1.30 

3.63 

± 

2.18 

-1.11 

+ 

0.09 

0.16 

+ 

0.01 

I Zw 191 

5.78 

± 

1.38 

9.27 

± 

1.76 

-1.28 

+ 

0.10 

0.57 

+ 

0.04 

IV Zw 93 

2.25 

± 

1.88 

6.89 

± 

2.02 

-1.21 

+ 

0.07 

0.22 

+ 

0.03 

Zw 2220 

6.74 

± 

2.37 

8.24 

± 

3.18 

-0.41 

+ 

0.12 

0.32 

+ 

0.06 

IV Zw 149 

1.49 

± 

3.11 

2.89 

± 

5.72 

0.02 + 0.21 

0.28 

+ 

0.01 

^'Gorrected 

for Galactic foregrc 

)und reddening 

1 I 2 OIII'). 

based or 

1 ISchlegel et alJ 

(Il998l') andlSchlaflv Sz Finkbeined 





= ln((HQ/H/3)/2.86), following the definition of ICalzetti et all 

il994) 
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log(0/H) + 12 log(0/H) + 12 log(0/H) + 12 


Fig. 2.— The SiIV_1400 (top) and CIV_1550 (bottom) indices measured from our HST 
spectra and compared to the KK04 (left). 03N2 (cent er), and N2 (right) metallicities 


determined from the optical data of iKong et al.l (j2002l ) . Solid lines illustrate the best 


linear fits, and Spearman’s p and p-values are given for each pair of diagnostics. 


metallicities. A positive relation between index value and metallicity is therefore slightly 
“flatter” under this definition than it would be for a metallicity-corrected continuum defi¬ 
nition (with the opposite holding through for a negative relation). For this reason we have 
measured the index values in our Starburst99 models from the non-normalized synthetic 
spectra, retaining the effects of a metallicity-dependent index continuum. These synthetic 
spectra are fully blanketed and are not adjusted for knowledge of the true continuum. Thus 
the same bias is present in both our observed and theoretical spectra. 


The evolution of the SiIV_1400 and CIV_1550 indices from our Starburst99 models as 
a function of metallicity are shown in Figures 3 and 4. Figure 3 shows models adopting 
an instantaneous burst star formation history at 1, 3, 5, and 7 Myr, while Figure 4 shows 
models adopting a continuous star formation history at an age of 5 Myr. In both cases 
the models are compared to this indices measured in this wor k, as shown in Figure 2, as 
well as the sub-sample of galaxies with measured indices from iLeitherer et al.l (120111 ) that 
are purely star-forming galaxies (i.e., no active galactic nucleus contribution). Our data 
assumes KK04 metallicities for this comparison; the choice of strong line diagnostic has 
only a minimal effect on the overall agreement between the data and the models. 
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Fig. 3.— SiIV_1400 (top) and CIV_1550 (bottom) indices measured from the SB99 models 
(black), simulating an instantaneous burst star formation history at (from left to right) 1, 
3, 5, and 7 Myr and adopting the M94 (circles), NOROT (triangles), and ROT (squares) 
stellar evolutionary tracks. The models are compared with the SiV_1400 and CIV_1550 


i ndice s measured from our HST data (blue) and the star-forming galaxies in lLeitherer et al. 
(1201 ih (red). All models assume a Kroupa IMF. 



Fig. 4.— As in Figure 4, but for models simulating a continuous burst star formation 
history at an age of 5 Myr. 
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Figures 3 and 4 show generally good agreement between the observations and models. 
For the SiIV_1400 index, the Starburst99 models generally underestimate the equivalent 
widths, particularly in the case of a continuous star formation history and the youngest ages 
of a single coeval stellar population (agreement improves at later ages). By contrast, in the 
case of the CIV_1550 index there is strong agreement between the models and observations 
for both the individual index values and the slope of the relation as a whole across all ages 
and star formation histories. The underestimation of the SiIV_1400 index by the models is 
ex pected; Starbur s t99 o nly simulates the stellar components of these indices, and as noted 


m 


Leitherer et al.l ([20111) interstellar components could account for as much as ~2/3 of the 


SiIV_1400 equivalent width. Therefore, the equivalent widths predicted by Starburst99 for 
the SiIV_1400 indices should be considered only as lower limits. 

It is also clear that the ROT models predict weaker equivalent widths f or both indices, 
which is explained by an o verall shift to a hotter population of stars (see iLevesaue et al. 


2 OI 2 I : iLeitherer et al.ll2014l i. Si iv (and to a lesser extent C iv) are not the dominant 
ionization states of Si and C, and a hotter population of ionizing stars will lead to weaker 
lines as a result (in contrast, features of dominant species suc h as N V and O vi are stronger 
in models that include stellar rotation: ILeitherer et al.ll2014l ). 


5. The /3 i 8 Extinction Diagnostic 


In optical spectra, the Ha/H/3 Balmer decrement is the most common diagnostic for 
quantifying extinction effects. By contrast, in the UV, the slope of the continuum is the 
most widely-used extinction diagnostic. The UV spectral slope, /3, is determined by fitting 
the continuum to the function f(A) oc A^, where /(A) is the flux density, to a specified 
wavelength range. ICalzetti et al.l (199J) consider both long-wavelength (1250-2800A) and 
short-wavelength (1250“1700A) definitions of /3, finding robust correlations with (where 
= In() , taking 2.86 as the theoretical Balmer decrement from Osterbrock 1989 1. 


‘B 


2.86 


Similarly, Calzetti ( 2001 1 defines ^ 26 , a ht spanning 1250-2600A, as a standardized measure 
of the UV slope and extinction diagnostic. 


Here we now consider /^ig, defined bv ICalzettil (|200ll ) as the continuum slope from 
I 25 O-I 8 OOA, as a potential additional diagnostic for shorter wavelengths (and, conse¬ 
quently, higher redshifts). Figure 5 shows the relationship between Pig and for our 
observed galaxies. In order to perform Pig fits to the UV continuua we first mask all of 
the Leitherer indices which fall within the diagnostic wavelength range. In the HST COS 
spectra, the geocoronal O i 1302 emission lines were also masked (the geocoronal Lya does 
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not need masking at the redshifts of our galaxies; see Figure 1). We calculated Ha/Hd 
ratios for each of our galax ies using the existing optical spectra from iKong et alJ (120021 ) , 
and define identically to lCalzetti et alJ (199J). To remove a ny potential conta mination 


from Galactic dust, we corrected both the raw line fluxes from 


Kong et alJ (120021') and ou r 


own UV spectra for Galactic reddening (ISchlegel et al.lll998l : ISchlafiv fc Finkbeineiil201ll ). 
Our values for /Jig and are given in Table 2, and our dat a show generally good agreement 
with the linear relationship determined in iGalzetti et alJ (|199J) for their definitions of a 
short-wavelength (3 and 


B- 



Fig. 5.— vs. Pis for our sample of obs erved galaxies. Values from are determined 
from the optical data of iKong et al.l (120(1211: values for dig ar e from our own HST spectra. 
The solid line is the best linear fit from lOalzetti et al.l ( 1994 ). The dashed lines represent 
the ~ 5(7 upper and lower envelopes of their data. 
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We note that ou r definition of /Jig spans a slightly different wavelength range than the 


Calzetti et al.l (199J) “short-wavelength” /J , covering 1250-1800A rather than the 1250- 


1700A. While the difference is slight this could account for some of the minor disagreements 
in Figure 5 between our data a nd the Calzetti relation (it is also worth noting that 16% of 
the original Calzetti et alJflQ^ data points also fall outside this relation). One galaxy, Haro 
25, does fall signihcantly above the rest of the sample; however, there are discrepancies 
in the literat u re reg arding the measured value of r from the optical spectrum. Data from 


Kehrig et al.l (120041 1 and the Hy/H/J decrement from Kong et al. (2002) both imply a 


higher reddening than the r determined here, which could partially improve this outlier’s 
agreement. 


/Jig holds particular value as an extinction diagnostic at high redshift. Spectral cov¬ 
erage in the UV is often limited to a narrower wavelength range than that spanned by 
the /J 26 diagnostic. At z > 3 this regime is redshifted into the near-IR; as a result, a UV 
extinction diagnostic over a shorter wavelength range would be advantageous for studying 
extinction in high-redshift galaxies. 


5.1. Theoretical Values of /Jig 

In Appendix A we provide tables of theoretical (dust-free) values for /Jig calculated 
from our Starburst99 models, to be used as references for determining extinction from 
both nearby and high-redshift galaxies using this dehnition of the UV continuum slope. 
Fully line-blanketed spectra were used to calculate the /Jig values, as opposed to using the 
true continuum. Our values are calculated for both instantaneous burst and continuous 
star-forming histories (at ages ranging from 0-20Myr in the instantaneous case and an 
equilibrium age of 5 Myr in the continuous case), all three stellar evolutionary tracks 
(M94, NOROT, and ROT), and all three IMFs simulated in our model grids (see Section 
3). 

The Starburst99 code is primarily based on stellar evolutionary synthesis, i.e., it is 
not a photoionization code. Therefore nebular physics is not considered and emission 
lines originating in HII regions and the general interstellar medium are not included in the 
SED. However, in the presence of hot ionizing stars continuous nebular emission can make a 
non-negligible contribution to the stellar emission and must be accounted for. Fortunately, 
nebular continuous emission is only weakly dependent on gas properties such as electron 
temperature (Tg) and electron density (rig), and detailed photoionization modeling is not 
required for estimating its contribution. Therefore this process is included in Starburst99. 
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Starburst99 in its default mode assumes all ionizing photons shortward of the Lyman 
break are absorbed in standard case B mode. We account for three continuous emission 
processes: bound-free -|- free-free emission of hydrogen; bound-free -|- free-free emission of 
neutral helium; and hydrogen two-phot on em i ssion . The absorption coefficients for these 
processes are taken from Tables 4-9 of I Alien (jl984l ): they are supplemented by values at 
longer wavelengths from iFerlandl ( 198fll ). Fixed values of Tg = 10,000 K and Ug = 10^ 
cm“^ are assumed (this assumption is non-critical for the parameter range of interest to 
this study). 



and collisional transitions to the 2p level occur, a transition from the 2s level to a virtual 
level takes place whose energy lies between those of the Is and 2s levels. This process is 
called two-photon emission since it produces two photons with a combined energy equal to 
that of the 2s to Is transition (i.e., Ly-a). The probability distribution of the process is 
symmetrical about half the energy level difference of 2s to Is transition, corresponding to 
a wavelength of 243lA. The peak emission of the two-photon spectrum is found to be at 
1621A after accounting for the weighting by the photon energy. 


Of the three nebular continuous emission processes included in Starburst99, the two- 
photon continuum by far dominates in the UV. The two-photon continuous absorption 
coefficients are an or der of magnitude large r than the bound-free -|- free-free coefficients 
between 1200-2000A. IPrake fc Ulrichl ( 198ll ) provide examples of the relative strengths of 
these processes. 


In order to assess the importance of the nebular contribution to the UV continuum, 
we computed two sets of models. One set is the Starburst99 default discussed above. This 
corresponds to case B in the radiation-bounded regime. Alternatively, we considered a 
second set having no nebular contribution. This could be a radiation-bounded nebula or a 
galaxy whose interstellar medium is porous due to the effects of stellar winds and super¬ 


novae which clear out cavities leading to the escape of ionizing photons (jBorthakur et al. 


20141 ). Comparison of the two model sets with and without nebular continuum suggests 


that omission of the nebular continuum results in a bluer spectrum (more negative /lig). 
The reason for this behavior is the redder spectral slope of the nebular compared with the 
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stellar continuum. Removal of the nebular component then leads to a bluer spectrum. 


Bouwens et al.l (|2010l l drew attention to the significance of the nebular continuum for 
the observed rest-frame UV spectra of galaxies at high redshift. The observed UV slopes 
are sometimes bluer than predicted by the most extreme theoretical models. Depending 
on dust reddening, metallicity, age, or initial mass function, the bluest slopes in the models 
have /3 > —2.7, whereas the observations in some galaxies suggest /3 < —3. Although 
observational uncertainties are lar ge and part of the d iscrepancies may be ascribed to larger 
than assumed photometric errors (iDunlop et al.ll2012l b th ere is tentative evi dence for bluer 
observed slopes when larger samples are considered (e.g.. I.Tiang et al.ll2ni3l ). Star-forming 
galaxies with redshifts approac hing the epoch of cosmic reionizat ion have been postulated 
to be leaking ionizing photons (jKuhlen fc Faucher-Giguere 2012 ). If so, one would expect 
a transition from the radiation-bounded to the density-bounded regime, and therefore a 
diminished contribution by the nebular continuum to the UV spectrum. 

Our two sets of /3is values given in Appendix A are thus intended to bracket the 
two extreme cases of 0 and 100 % escape frac tion (/esc) of ionizing photons. They are 
complementary to the models of llnouel (20 111) who provided a grid of nebular models, 
including continuum models, using an older version of Starburst99 when evolution models 
with rotation had not yet been available. 


The tables in Appendix A indicate several trends of the behavior of Pig. There is an 
overall evolution towards redder slopes with age, which is simply the result of on average 
lower stellar temperatures with older age. As an example, see Figure 6. The non-monotonic 
behavior around 4 to 6 Myr is due to the first appearance of luminous post-main-sequence 
stars whose somewhat cooler temperature leads to a redder UV continuum. These stars 
do not subsequently evolve into red supergiants but rather turn into hot Wolf-Rayet stars. 
Overall, the variation of /3i8 with age is rather moderate, which of course is the reason for 
the suitability of the UV slope as a dust indicator. 


The three IMF choices considered here result in only minor differences in The 
Kroupa IMF has parameters {aiow = 1-3, ahigh = 2.3). An IMF enriched in massive stars 
relative to the Kroupa IMF (1.3, 1.3) has a UV slope which is bluer by ~ 0.2 than an IMF 
deficient of massive stars (1.3, 3.3). More extreme variations of /3i8 are expected in the 
presence/absence of very massive stars with masses higher than, e.g., 150 Mq. These stars 
are not covered by the evolutionary tracks included in Starburst99. 

The two extreme cases of /esc = 0 and 100 (all ionizing photons absorbed or escaping, 
respectively) have a rather significant effect on /3i8 for ages < 10 Myr. For ages and IMF 
values favoring the production of ionizing radiation, the UV slopes with /esc = 100 are 
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Fig. 6.— /?i8 measured from the SB99 models, simulating an instantaneous burst star 
formation history with ROT evolutionary tracks and a Kroupa IMF. Metallicities are Z = 
0.002 (red) and Z = 0.014 (blue). The two extreme cases of /esc = 0 (left) and 100 (right) 
are illustrated. 

bluer by ~ 0.5. This behavior is easily understandable from the previous discussion of the 
nebular continuum: the nebular continuum is dominated by the two-photon process whose 
slope is redder than that of the stellar continuum. 

Metallicity has a relatively minor impact on /3i8. While individual metal-poor stars 
have bluer colors than metal-rich stars, stellar evolution at high metallicity favors the 
presence of hot Wolf-Rayet stars, which compensates for the color effect in individual 
stars. As a result, the /3i8 values at various metallicities are not very different. We remind 
the reader that this argument only applies to models for single-star evolution. If significant 
numbers of Wolf-Rayet stars at low metallicity are produced in binaries, /3i8 may become 
bluer as well. 

Models with rotation generally result in more negative /3i8 values than non-rotating 
models. The effect is less pronounced for /esc = 0 because rotation leads to hotter stars 
which are bluer on the one hand but are also responsible for a stronger nebular continuum 
on the other. 
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6. Discussion and futnre work 


We have demonstrated that there is a statistically robust positive linear correlation 
with metallicity for the UV stellar wind indices SiIV_1400 and CIV_1550. This correlation 
is supported both by multi-wavelength spectroscopic observations of star-forming galaxies 
and by the predictions of the Starburst99 stellar population synthesis models. This suggests 
that parameterizations of UV metallicity diagnostics based on synthetic galaxy spectra will 
offer a well-calibrated suite of new global metallicity diagnostics that can be applied to 
spectra of star-forming galaxies at high redshifts. 

In addition, the clear correlation between the UV and optical metallicity diagnostics 
has important implications for future studies of star-forming galaxy metallicities across a 
broad range of redshifts. While a larger and higher-S/N multi-wavelength observational 
sample is required to fully sample the parameter space of ISM properties and diagnos¬ 
tics in star-forming galaxies (including star formation rate, young stellar population age, 
stellar mass, and extinction), our results demonstrate that it will be possible to compare 
metallicities determ ined from both UV and optical spectra through conversions similar to 
those calculated bv iKewlev &: EllisonI ()2008l i for optical diagnostics alone. This will offer 
the first quantitative means of unifying metallicity measurements for star-forming galaxies 
across multiple wavelengths, and consequently across a much broader range of redshifts. 

Similarly, combining the ionizing SED outputs of the full grid of Starburst99 rnodels 

(including those acc ommodating the effects of stellar rotation, e.g. Levesque et al. 20121: _ 

Leitherer et al.ll2014) wit h a photoionizat ion code such as MA PPINGS (e.g. Sutherland Sz Donita 


1993 


Groves et al 


20041) or CLOUDY ( Ferland et al.lll998l) would allow us to determine 


theoretical parameterizations for indices with strong or dominant contributions from inter¬ 
stellar lines and directly compare these to predicted fluxes for optical emission lines. This 
includes the interstellar components of the SiIV_1400 and CIV_1550 indices as well as other 
far-UV diagnostic indices that are primarily interstellar such as CII_1330 and A1II_1670. 

We have also determined theoretical values of Pis that can be used for determining 
extinction in rest-frame UV spectra of star-forming galaxies, assuming the two extremes of 
0 % and 100% /esc to facilitate their use with observations of star-for ming galaxies in both 
the nearby and hig h-redshift universe. Following previous work by ICalzetti et al.l (|l994l ) 
and ICalzettil (j200lh , Pis shows a similar correlation with optical extinction diagnostics and 
is particularly valuable for high-redshift observations. In the future, a combination of stellar 
population synthesis and photoionization models would also allow for a valuable theoretical 
comparison between the extinction diagnostics Pis and r^, allowing us to directly probe 
the effects of other galaxy properties (such as stellar population age, ionization parameter. 































and electron density) on these diagnostics and improving their utility at high redshift. 


This paper was based on observations made with the NASA/ESA Hubble Space Tele¬ 
scope associated with program 13481; support for this program was provided by NASA 
through a grant from the Space Telescope Science Institute. EML is additionally supported 
by NASA through Hubble Eellowship grant number HST-HE-51324.01-A from the Space 
Telescope Science Institute. The Space Telescope Science Institute is operated by the As¬ 
sociation of Universities for Research in Astronomy, Incorporated, under NASA contract 
NAS5-26555. 
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A. /3i 8 in Star-Forming Galaxy Models 


This preprint was prepared with the A AS lAERK macros v5.2. 



Table Al. /Jig in Starburst99 Models with M94 Evolutionary Tracks, /esc = 0 


Age _IMF = 1.3,1.3_IMF = 1.3,2.3_IMF = 1.3,3.3 


(Myr) 

0.001 

0.004 

0.008 

0.020 

0.040 

0.001 

0.004 

0.008 

0.020 

0.040 

0.001 

0.004 

0.008 

0.020 

0.040 

0.0 

-2.3582 

-2.4315 

-2.4681 

-2.5303 

-2.6174 

-2.4936 

-2.5563 

-2.5763 

-2.5999 

-2.6655 

-2.6285 

-2.6516 

-2.6432 

-2.5977 

-2.6261 

2.0 

-2.3579 

-2.4989 

-2.5528 

-2.5164 

-2.4756 

-2.4870 

-2.5813 

-2.6138 

-2.5888 

-2.5010 

-2.6258 

-2.6460 

-2.6365 

-2.6114 

-2.5195 

4.0 

-1.9900 

-2.3550 

-2.1951 

-2.0406 

-1.7666 

-2.2533 

-2.4907 

-2.3632 

-2.1892 

-1.9992 

-2.5178 

-2.5926 

-2.5014 

-2.3348 

-2.2026 

6.0 

-2.3226 

-2.1695 

-2.1999 

-2.1696 

-2.4203 

-2.4730 

-2.3546 

-2.3646 

-2.3147 

-2.4869 

-2.5839 

-2.4851 

-2.4670 

-2.3974 

-2.4741 

8.0 

-2.4624 

-2.4085 

-2.7041 

-2.6920 

-2.6077 

-2.5612 

-2.5170 

-2.7184 

-2.6844 

-2.5892 

-2.6103 

-2.5528 

-2.6493 

-2.5932 

-2.4937 

10.0 

-2.4620 

-2.4402 

-2.7933 

-2.7155 

-2.5876 

-2.5528 

-2.5290 

-2.7694 

-2.6807 

-2.5460 

-2.5864 

-2.5392 

-2.6561 

-2.5602 

-2.4281 

12.0 

-2.4640 

-2.4698 

-2.7935 

-2.6916 

-2.5423 

-2.5471 

-2.5368 

-2.7524 

-2.6398 

-2.4853 

-2.5651 

-2.5193 

-2.6203 

-2.5031 

-2.3540 

14.0 

-2.4446 

-2.4952 

-2.7643 

-2.6375 

-2.4646 

-2.5227 

-2.5396 

-2.7124 

-2.5776 

-2.4009 

-2.5319 

-2.4963 

-2.5703 

-2.4351 

-2.2652 

16.0 

-2.4502 

-2.4944 

-2.6345 

-2.3670 

-2.3853 

-2.5173 

-2.5238 

-2.6003 

-2.3579 

-2.3176 

-2.5114 

-2.4655 

-2.4767 

-2.2663 

-2.1803 

18.0 

-2.4352 

-2.3870 

-2.3452 

-2.2390 

-2.3084 

-2.4940 

-2.4275 

-2.3702 

-2.2421 

-2.2375 

-2.4795 

-2.3845 

-2.3117 

-2.1665 

-2.0993 

20.0 

-2.4114 

-2.3583 

-2.3309 

-2.2279 

-2.2330 

-2.4648 

-2.3915 

-2.3434 

-2.2154 

-2.1597 

-2.4447 

-2.3428 

-2.2744 

-2.1254 

-2.0207 

5.0 (cont) 

-2.2595 

-2.3961 

-2.3491 

-2.3757 

-2.2533 

-2.4101 

-2.5119 

-2.4617 

-2.4431 

-2.3015 

-2.5851 

-2.6115 

-2.5573 

-2.4966 

-2.3722 


KS 

Table A2. /Jig in Starburst99 Models with M94 Evolutionary Tracks, /esc = 100 i 


Age _IMF = 1.3,1.3_IMF = 1.3,2.3_IMF = 1.3,3.3 


(Myr) 

0.001 

0.004 

0.008 

0.020 

0.040 

0.001 

0.004 

0.008 

0.020 

0.040 

0.001 

0.004 

0.008 

0.020 

0.040 

0.0 

-3.0116 

-3.0182 

-3.0274 

-3.0340 

-3.0690 

-2.9980 

-3.0020 

-2.9979 

-2.9767 

-2.9989 

-2.8785 

-2.8718 

-2.8553 

-2.7920 

-2.7996 

2.0 

-3.0053 

-2.9722 

-2.9671 

-2.7910 

-2.6393 

-2.9944 

-2.9667 

-2.9464 

-2.8199 

-2.6421 

-2.8914 

-2.8603 

-2.8247 

-2.7521 

-2.6135 

4.0 

-2.1098 

-2.4927 

-2.3095 

-2.1159 

-1.8247 

-2.4001 

-2.6240 

-2.4727 

-2.2560 

-2.0443 

-2.6481 

-2.6930 

-2.5851 

-2.3856 

-2.2338 

6.0 

-2.4094 

-2.2194 

-2.2319 

-2.1838 

-2.4660 

-2.5616 

-2.4063 

-2.3972 

-2.3295 

-2.5198 

-2.6549 

-2.5274 

-2.4937 

-2.4103 

-2.4945 

8.0 

-2.5034 

-2.4293 

-2.7203 

-2.6993 

-2.6113 

-2.6016 

-2.5370 

-2.7328 

-2.6910 

-2.5925 

-2.6429 

-2.5686 

-2.6600 

-2.5983 

-2.4964 

10.0 

-2.4773 

-2.4479 

-2.8002 

-2.7189 

-2.5893 

-2.5679 

-2.5365 

-2.7755 

-2.6838 

-2.5476 

-2.5988 

-2.5453 

-2.6606 

-2.5626 

-2.4293 

12.0 

-2.4708 

-2.4738 

-2.7972 

-2.6934 

-2.5432 

-2.5538 

-2.5406 

-2.7556 

-2.6415 

-2.4861 

-2.5707 

-2.5223 

-2.6227 

-2.5043 

-2.3546 

14.0 

-2.4483 

-2.4978 

-2.7666 

-2.6387 

-2.4651 

-2.5263 

-2.5419 

-2.7144 

-2.5786 

-2.4014 

-2.5349 

-2.4981 

-2.5718 

-2.4359 

-2.2656 

16.0 

-2.4525 

-2.4961 

-2.6360 

-2.3677 

-2.3857 

-2.5196 

-2.5253 

-2.6016 

-2.3584 

-2.3179 

-2.5132 

-2.4667 

-2.4777 

-2.2667 

-2.1805 

18.0 

-2.4368 

-2.3880 

-2.3459 

-2.2393 

-2.3086 

-2.4956 

-2.4285 

-2.3709 

-2.2425 

-2.2378 

-2.4808 

-2.3853 

-2.3123 

-2.1668 

-2.0995 

20.0 

-2.4125 

-2.3590 

-2.3314 

-2.2282 

-2.2332 

-2.4659 

-2.3921 

-2.3439 

-2.2156 

-2.1598 

-2.4457 

-2.3434 

-2.2748 

-2.1256 

-2.0208 

5.0 (cont) 

-2.6150 

-2.6999 

-2.6211 

-2.6008 

-2.4184 

-2.7258 

-2.7688 

-2.6842 

-2.6152 

-2.4231 

-2.7844 

-2.7695 

-2.6938 

-2.5996 

-2.4448 
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Table A3. /?i8 in Starburst99 Models with NOROT Evolutionary Tracks, /esc = 0 


Age 

(Myr) 

IMF = 

1.3,1.3 

IMF = 

1.3,2.3 

IMF = 

1.3,3.3 

0.002 

0.014 

0.002 

0.014 

0.002 

0.014 

0.0 

-2.4064 

-2.5104 

-2.5245 

-2.5905 

-2.6395 

-2.6198 

2.0 

-2.4144 

-2.5906 

-2.5271 

-2.6468 

-2.6431 

-2.6423 

4.0 

-2.0675 

-2.0281 

-2.2988 

-2.1726 

-2.5236 

-2.3276 

6.0 

-2.2837 

-2.1980 

-2.4284 

-2.3378 

-2.5363 

-2.4155 

8.0 

-2.4082 

-2.6100 

-2.5132 

-2.6303 

-2.5630 

-2.5599 

10.0 

-2.4129 

-2.6489 

-2.5098 

-2.6357 

-2.5378 

-2.5282 

12.0 

-2.4061 

-2.6926 

-2.4928 

-2.6423 

-2.5039 

-2.4993 

14.0 

-2.3747 

-2.6771 

-2.4548 

-2.6068 

-2.4569 

-2.4466 

16.0 

-2.3310 

-2.6169 

-2.4080 

-2.5405 

-2.4058 

-2.3770 

18.0 

-2.3233 

-2.5536 

-2.3907 

-2.4739 

-2.3765 

-2.3097 

20.0 

-2.3581 

-2.4955 

-2.4071 

-2.4131 

-2.3723 

-2.2489 

5.0 (cont) 

-2.3203 

-2.3289 

-2.4545 

-2.4249 

-2.6024 

-2.5026 


Table A4. /3i8 in Starburst99 Models with NOROT Evolutionary Tracks, /esc = 100 


Age 

(Myr) 

IMF = 

1.3,1.3 

IMF = 

1.3,2.3 

IMF = 

1.3,3.3 

0.002 

0.014 

0.002 

0.014 

0.002 

0.014 

0.0 

-3.0202 

-3.0367 

-3.0071 

-2.9926 

-2.8862 

-2.8306 

2.0 

-3.0074 

-2.9072 

-3.0017 

-2.9155 

-2.9001 

-2.8042 

4.0 

-2.1767 

-2.1083 

-2.4284 

-2.2469 

-2.6385 

-2.3862 

6.0 

-2.3406 

-2.2177 

-2.4872 

-2.3579 

-2.5854 

-2.4325 

8.0 

-2.4278 

-2.6183 

-2.5332 

-2.6378 

-2.5798 

-2.5656 

10.0 

-2.4208 

-2.6528 

-2.5176 

-2.6391 

-2.5444 

-2.5308 

12.0 

-2.4100 

-2.6948 

-2.4966 

-2.6442 

-2.5071 

-2.5007 

14.0 

-2.3771 

-2.6785 

-2.4571 

-2.6080 

-2.4589 

-2.4475 

16.0 

-2.3326 

-2.6177 

-2.4095 

-2.5413 

-2.4071 

-2.3775 

18.0 

-2.3245 

-2.5542 

-2.3918 

-2.4744 

-2.3774 

-2.3101 

20.0 

-2.3590 

-2.4959 

-2.4081 

-2.4135 

-2.3730 

-2.2492 

5.0 (cont) 

-2.6522 

-2.5524 

-2.7447 

-2.6072 

-2.7865 

-2.6166 
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Table A5. /3i8 in Starburst99 Models with ROT Evolutionary Tracks, /esc = 0 


Age 

(Myr) 

IMF = 

1.3,1.3 

IMF = 

1.3,2.3 

IMF = 

1.3,3.3 

0.002 

0.014 

0.002 

0.014 

0.002 

0.014 

0.0 

-2.4215 

-2.5341 

-2.5448 

-2.6020 

-2.6486 

-2.5968 

2.0 

-2.4508 

-2.5753 

-2.5539 

-2.6382 

-2.6527 

-2.6316 

4.0 

-2.2227 

-2.2496 

-2.3848 

-2.3583 

-2.5619 

-2.4609 

6.0 

-2.1270 

-2.0764 

-2.3188 

-2.2163 

-2.4825 

-2.3388 

8.0 

-2.1814 

-2.3088 

-2.3377 

-2.3892 

-2.4586 

-2.4229 

10.0 

-2.3127 

-2.6991 

-2.4309 

-2.6865 

-2.4946 

-2.5940 

12.0 

-2.4371 

-2.5492 

-2.5180 

-2.5628 

-2.5330 

-2.4866 

14.0 

-2.4955 

-2.6723 

-2.5533 

-2.6280 

-2.5402 

-2.4979 

16.0 

-2.4309 

-2.6436 

-2.4918 

-2.5861 

-2.4814 

-2.4441 

18.0 

-2.3848 

-2.5977 

-2.4446 

-2.5323 

-2.4325 

-2.3838 

20.0 

-2.3579 

-2.5487 

-2.4131 

-2.4780 

-2.3949 

-2.3260 

5.0 (cont) 

-2.3655 

-2.3248 

-2.4828 

-2.4522 

-2.6107 

-2.5318 


Table A6. j3i% in Starburst99 Models with ROT Evolutionary Tracks, /esc = 100 


Age 

(Myr) 

IMF = 

1.3,1.3 

IMF = 

1.3,2.3 

IMF = 

1.3,3.3 

0.002 

0.014 

0.002 

0.014 

0.002 

0.014 

0.0 

-3.0191 

-3.0324 

-3.0043 

-2.9739 

-2.8775 

-2.7874 

2.0 

-3.0005 

-3.0328 

-2.9970 

-2.9822 

-2.8940 

-2.8166 

4.0 

-2.3642 

-2.5496 

-2.5529 

-2.5959 

-2.7096 

-2.6085 

6.0 

-2.2135 

-2.1978 

-2.4089 

-2.3166 

-2.5574 

-2.4081 

8.0 

-2.2183 

-2.3523 

-2.3755 

-2.4243 

-2.4906 

-2.4473 

10.0 

-2.3295 

-2.7125 

-2.4474 

-2.6972 

-2.5084 

-2.6015 

12.0 

-2.4456 

-2.5527 

-2.5263 

-2.5660 

-2.5398 

-2.4891 

14.0 

-2.5003 

-2.6745 

-2.5579 

-2.6300 

-2.5439 

-2.4994 

16.0 

-2.4337 

-2.6449 

-2.4945 

-2.5873 

-2.4837 

-2.4449 

18.0 

-2.3866 

-2.5986 

-2.4464 

-2.5331 

-2.4340 

-2.3843 

20.0 

-2.3593 

-2.5493 

-2.4145 

-2.4786 

-2.3960 

-2.3265 

5.0 (cont) 

-2.7260 

-2.6853 

-2.7964 

-2.7344 

-2.8079 

-2.6949 




